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Abstract 

We revisit a model in which the ionization energy of a metal particle is associated with the work 
done by the image charge force in moving the electron from infinity to a small cut-off distance just 
outside the surface. We show that this model can be compactly, and productively, employed to 
study the size dependence of electron removal energies over the range encompassing bulk surfaces, 
finite clusters, and individual atoms. It accounts in a straightforward manner for the empiri- 
cally known correlation between the atomic ionization potential (IP) and the metal work function 
(WF), IP/WF~2. We formulate simple expressions for the model parameters, requiring only a 
single property (the atomic polarizability or the nearest neighbor distance) as input. Without any 
additional adjustable parameters, the model yields both the IP and the WF within ^^10% for all 
metallic elements, simulates the concentration dependence of the WF of regular binary bulk alloys, 
and matches the size evolution of the IP of finite metal clusters for a large fraction of the experi- 
mental data. The parametrization takes advantage of a remarkably constant numerical correlation 
between the nearest-neighbor distance in a crystal, the cube root of the atomic polarizability, and 
the image force cutoff length. The paper also includes an analytical derivation of the relation of 
the outer radius of a cluster of close-packed spheres to its geometric structure. 

PACS numbers: 79.60.Jv, 78.67.-n, 61.46. +w 



I. INTRODUCTION 



While the good agreement between theoretical and experimental atomic ionization po- 
tentials (IP) is a major triumph for quantum mechanics, it is prohibitively more difficult 
to rigourously solve the polyatomic quantum problem, not to mention extrapolation to an 
infinite bulk metal. The IP for an atom is a well understood quantity. The same cannot, 
however, be said about the work function (WF) for a metal or even for a finite cluster. On 
the other hand, it is an experimentally realistic task to produce clusters ranging in size from 
two atoms up to tens of thousands of atoms and to measure the size dependence of the elec- 
tron removal energy. The clusters can be made big enough that they can be considered close 
to bulk metals, hence the evolution from the atomic IP to the metal WF can be mapped 
out. Experiments over such a wide range have been performed, e.g., for sodiumB. For each 
element in the periodic table one would therefore expect that there exists a function which 
can predict the electron removal energy for a particle of arbitrary size. 

The exact derivation of such a scaling law is a daunting task. However, the available 
experimental data on the IP, the WF, and on clusters reveal some characteristic features. 
For example, it has been noted a long time ago that the IP and the WF of metallic elements 
are approximately correlated to each other via the factoiil 

I/W ^ 2. (1) 

To give another example, for many metal clusters of intermediate sizes the electron removal 
energy has been found to scale as 

0^W^ + 7^, (2) 

where R is the radius of the cluster and 7 is a constant factoi^S. Henceforth, / will denote 
the atomic first ionization potential, W will denote the polycrystalline bulk surface work 
function, and for a finite metal particle the term "electron removal energy" will be employed 
(denoted by 0). 

A significant number of first-principles calculations have been performed for work func- 
tions of metallic systemsi'iJi. In addition, simple models based on semi-empirical approaches 
combined with classical electrostatics have been rather successful at reproducing WF trends 
and values!. This suggests that some features of the desired scaling law may be found by 
employing such a model. In this paper we demonstrate that by combining the image-charge 
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potential function for a finite particle with just a single material-dependent scale parameter 
(the atomic polarizability a or the crystalline nearest-neighbor distance r„„) one can obtain 
an interpolation formula covering the full size range from the atom through the cluster to 
the bulk. This formula estimates both the IP and the WF within ~ 10% for all metallic 
elements in the periodic table, yields values in reasonable agreement with experiment for 
many intermediate sized clusters, and provides a natural justification for the aforementioned 
IP/WF^ 2 ratio. 

The plan of the paper is as follows. In Sec. || we consider the image-charge expression 
which describes the removal of an electron from an isolated sphere. By focusing on the limits 
of a sphere of atomic radius and one of infinite radius, we show that the relation (JlD follows as 
naturally from the image-charge consideration as does Eq. (0). In Sec. |Tl| we demonstrate 
some striking parallels between the variation of atomic sizes and image-force cutoff distances 
across the periodic table and use these observations to formulate compact expressions for 
estimating both WF and IP. An extension to the case of binary alloys is discussed in Sec. 

Finally, we show that by using an appropriately interpolated expression for the cluster 
radius a good description of the size evolution of cluster ionization potentials can be obtained 
(Sec. |V]). Some rigorous formulas on close-packed cluster radii are derived in the Appendix. 



II. THE IMAGE POTENTIAL MODEL 



It is stated in the literature that the earliest attempt to explain the work function of a 
metallic surface using classical electrostatics was due to Deby(Jli. He proposed that it is 
equal to the energy required to pull an electron out to infinity against its image charge. 
Since the image force diverges at the surface, Schottkylll suggested that one may be able to 
define a microscopic cutoff distance d at which the image force starts to act. Despite the 
simplicity of this modeS, there have been numerous attempts to estimate the parameter d 
in order to fit the experimental values of the WFii'ilii00. Although the particular choices 
of d were supported only by plausibility arguments, they were frequently able to offer rather 
nice agreement with the experimental data. 

In a similar spirit, let us now consider the image- force expression for the energy required 
to remove an electron from an isolated finite metal particle, modelled as a conducting sphere 
of radius R. The particle is assumed initially neutral, i.e., after the removal of the electron it 
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acquires an unit positive charge. A calculation of the work required to move the electron from 
a distance d outside the metal surface to infinity against its image charge is a straightforward 
exerciselll: 

_ ( l+A{d/R) + (S{d/RY + 2{dlRf \ 
^^^^ ~ Ad \ {l + d/2R){l + d/RY ) ■ 
The cutoff parameter d is assumed to be a material-dependent constant. The first factor on 
the right-hand side represents the bulk {R — oo) work function: W = /{Ad). This formula 
has been applied to yttrium and lanthanide clusters in Ref. |T^. However, the authors did 



not extend it to the bulk or the atomic limit; they used the equation as an extrapolation 
formula for small clusters with the WF as a given boundary condition. 
It is convenient to rewrite the above equation in the following form: 

^ = ndiR), (4) 

where the dimensionless scaling function ri{d/R) is defined by the second factor in Eq. @ 
and plotted in Fig. |l|. If Eq. were applied all the way down to the atomic limit, it would 
give an estimate of the ratio IP/WF as the value of ri{d/Rat), where Rat is a quantity charac- 
terizing the atomic size. On the other hand, we know from numerous investigationsi'ii'Ei'Blli 
that the cutoff parameter d is, sensibly enough, also of the same order of magnitude. In the 
atomic limit, therefore, the ratio d/Rat should be on the order of unity. In other words, if 
one assumes that d ~ Rat, then, independent of the exact expression for either parameter, 
the scaling function predicts that 

I/W ~ 7]{1) ^ 2. (5) 

This is a new explanation of the well known empirical result mentioned in the Introduction, 
Eq. (|l|). In the next section we suggest some specific parametrizations of the variables d 
and Rat and show that these can give an even more accurate value of the IP/WF ratio. 

In the large R limit, the scaling function can be expanded to first order in d/R. The 
result is: 

,(«).H'.?^H-0((|)^). (6) 

This is the well known finite size correction for the ionization potential of metallic 
clusterli'Sffl. This scaling law has been experimentally verified for many metal 
clusterlS§i. Although there is still some controversy whether the 3/8 factor is sufficiently 
rigorousii'0'0@ii, Eq. (^) does fit the experimental data relatively well. 



It appears, therefore, that Eq. (j^) offers a consistent estimate for the scahng of electron 
removal energy from the atomic IP to finite particle sizes to the bulk WF. It is interesting 
to ask whether some simple parametrizations for the image force cutoff distance and atomic 
and cluster radii may be proposed so as to enable more quantitative applications of Eq. 
to experimental data. This is the subject of the sections that follow. 



III. LENGTH SCALE PARAMETERS 



For guidance with length scaling, let us begin by plotting the "experimental" image force 
cutoff distances as defined by equating e^/ {4d) to the experimentally measured bulk surface 
polycrystalline work functions^. These values are shown in Fig. H(a). The same figure 
displays two quantities reflecting the size of individual atoms: the cube root of the atomic 
polarizabilityill (r^) and one-half of the nearest neighbor distance in a crystalil (r„„/2). 
Although d shows strong correlation with both the polarizability radius and the nearest 
neighbor distance, it's consistently lower than either one. However, by forming the inverse 
sum of the two, a quantity emerges which follows the experimental values of d beautifully. 
This parametrization has the form 

1/d ^ l/r„ + 2/rnn (7) 

and agrees with the empirical value of the cutoff parameter to within ~ 10% for most 
metallic elements, see Fig. 0(b). Put another way, the WF for most metallic elements can 
be estimated to better than 10% by using this parametrization of d. 

The coexistence of an atomic part and a bulk part in a work function model resembles 
qualitatively the elegant argument in Ref. in which the WF value is calculated as arising 



from the IP of a neutral atom reduced by the work done by the image force in bringing the 
resulting ion back to the crystal surface. 

As a matter of fact. Fig. |^(a) reveals that all the size parameters turn out to have 
essentially the same trend across the periodic table, differing only in overall magnitude. 
We are not aware of a quantitative theory explaining this observation for the crystalline 
state, but, roughly speaking, one does expect bigger atoms to have higher polarizabilities, 
larger lattice spacings in crystal form, and lower IP together with lower WF [via Eqs. 
and therefore greater d values. In fact, remarkably linear correlations between the 
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polarizability radius and such properties of individual atoms as the radius of the outermost 
orbital, the ionization potential, and the electronegativity have been extensively discussed 
in the literature (see Refs. pO| , pT| , p2| and references therein). This strongly suggests that just 
a single size parameter may be sufficient to yield an estimate for the WF. 

Indeed, as shown in Fig. ^, it turns out that the ratio r„„/(2rQ) is very close to 0.65 
for all metallic elements^ (the deviation from this average is within ~ 20%), and this 
correlation can be substituted into Eq. (0). In the end, therefore, one can use, e.g., only the 
polarizability radius in order to predict the approximate value of the bulk surface WF. 

As shown in Fig. |, a plot of 

= (8) 
i.e., Eq. (|^) with r„„/2 = 0.65rQ, [cf. Fig. 0(c)] confirms that this single-parameter formula 
gives WF estimates within ~ 10% of the experimental values. As mentioned at the beginning 
of Sec. H, a variety of expressions for the quantity d have been proposed earlier, but 
the one given here shares the virtues of needing only a single piece of information (the 
atomic polarizability or the nearest neighbor distance), being free of additional adjustable 
parameters, and producing rather accurate results. 

Returning to the subject of the IP/WF ratio discussed in the previous section, Eq. (^, 
we notice that even more accurate agreement with the experimental average is obtained if 
we employ 

d/Rat ^ r„„/(2r,) ^ 0.65. (9) 

Indeed, this is a reasonable representation of the ratio of a characteristic crystalline length 
to an atomic length. This refinement gives 

I/W ^ r/(0.65) = 1.8, (10) 

see Fig. ||. 

With estimates of both the work function WF, and the ratio IP/WF ^ 1.8, the atomic 
first ionization potential can be calculated. Fig. ^ shows the calculated IP plotted vs. the 
exact known IP0. The accuracy is again ± 10% for most elements^. 

Eqs. (^ and (^) imply the following expression for the size parameter Rat in the atomic 
limit: 

Rat ~ .""^ ,^ ~ 0.61r„ (11) 



where in the last part we again used the average experimental ratio between r„„ and r^. 

The discussion in this section shows that a consistent and economical description of 
electron removal energies both in the bulk and the atomic limits can be achieved with the 
help of only a single atomic parameter. In the next section we will consider the region in 
between, i.e., the electron removal energies for clusters of a finite number of atoms. 



IV. ESTIMATE OF ALLOY WORK FUNCTIONS 

The regularities discussed above have been deduced for pure metals. It is interesting 
to inquire what they might imply for mixed systems such as binary alloys. Obviously, it is 
beyond the capacity of the present approach to achieve a microscopic understanding of mixed 
systems and to incorporate the possible variations in structure and miscibility . Nevertheless, 
it makes sense to consider a simple scaling argument for continuous bulk solutions. 

The preceding section demonstrated that the WF can be rather accurately parametrized 
as being inversely proportional either to the polarizability radius or to the nearest-neighbor 
distance. In an alloy AxBi_x the effective values of these length parameters will change. 
In fact, there are two natural approximations which can be made. One is that the lattice 
parameter will change linearly with concentration. This is known as Vegard's lawHS Within 
the framework of the present model, the alloy WF would then behave as the inverse average 
of the constituents' WF's: l/W{x) ^ x/Wa + (1 - x)/Wb- 

An alternative and possibly more physical assumption is that it is the effective polariz- 
ability that should be taken as varying linearly with concentration. Indeed, the polarizability 
of a small volume element of the sample is approximately the weighted average of the po- 
larizabilities of the components. In this case, with W ~ a^^^'^, we anticipate 

[W{x)f ^Wl^ 'W% ^ ' 

Note that this relation would also be obtained if we assumed that it was the molar 
volumes, rather than lattice parameters, that would need to be averaged upon alloying 
(Zen's lawii). 

Well-controlled measurements of W{x), like rigorous calculations, are not many in num- 
ber: traces of surface segregation or contamination can obliterate the basic effect. Fig. |^ 
illustrates the prediction of Eq. (|T2p for binary alloys of silver and gold, carefully charac- 



terized in Ref. ^ The scaling curve reproduces the trend of the data and even the correct 
curvature of the Ag-Au hne. It is noteworthy that predicting the sign of the deviation of 
W{x) from hnearity appears to be a difficult challenge for microscopic theory. The one 
calculation which focused on such a deviationil applied a self-consistent two-band charge- 
transfer model to Ag^Aui^x and was able to yield the correct sign only when level shifts due 
to alloying were taken into account. The absolute magnitude of the calculated deviation was 
approximately one-third of the experimental value, and was in fact very close to the curve 
in Fig. H produced by Eq. ([121) . 



For cases when the constituents' work functions differ by relatively small amounts, such 
as alkah alloys§'§, (Eq. ([12|) predicts an essentially linear interpolation, in agreement with 
the data. 

It is seen, therefore, that scaling estimates can be usefully applied to the work functions 
of smoothly evolving alloys. Of course, cases when the WF is not a monotonic function of 
concentrationSill call for more involved structural and theoretical analysis. 

V. ELECTRON REMOVAL ENERGIES FOR CLUSTERS AND HARD SPHERE 
PACKING 

Having shown that the scaling law, Eq. (^, offers good guidance for both IP and WF, we 
now address the intermediate regime of finite clusters. Therefore a workable parametrization 
of cluster radius, Rcu is needed for use in this equation. 

A commonly used method to estimate the particle radius is to scale it according to the 
Wigner-Seitz radius: Rci{N) = r^N^^^, where is the number of atoms. This expression, 
commonly augmented with an electron spill-out term, has been used to describe many cluster 
properties. However, if this formula is used in Eq. (|^) together with the expressions for d 
obtained above, a close match with the experimental data is not obtained. Moreover, it is 
desirable to retain a form that would depend only on A^ and on the parameters Tq, and/or 
Vnn, as in the atomic and bulk limits described above. 

A complementary definition of cluster radius has been proposed in Ref. |1T7| , ^ . The 



cluster was modelled as a group of hard spheres packed so as to minimize the surface area. 
The radius of a sphere circumscribing the close-packed cluster was then substituted into Eq. 
(0). By packing steel balls in a rubber envelope, the authors found empirically that this 
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circumscribing sphere radius could be expressed as 

RciiN) ^ l.SroN'^^ (13) 

for clusters N 7, where tq is the radius of each individual hard sphere. For smaller 
clusters, the exact circumscribing sphere radius can be derived from simple geometry^. This 
definition gave reasonable answers for several metal clusters, although additional adjustable 
parameters had to be employed. We found that an expression of the type ([T3|) can be derived 
analytically for various cluster packing arrangement. The precise numerical coefficient varies 
with the packing structure, but in all cases is not far from 1.3. The analysis is described in 
the Appendix. 

Based on these considerations, we formulate the following interpolation for the parameter 
R to be used in Eq. @ for the calculation of clusters' 4>{R): 

^ RUN) 

RUN) + r„„/2- ^ ' 

Here Rci{N) is the outer, or circumscribing (as described above) radius for a cluster with 
hard spheres, each of a radius equal to the polarizability radius ro = r^. As before, r„„ 

is the nearest neighbor distance of the bulk metal . For A^ ^ 1 this expression turns into 

rM, while in the atomic limit it is the same as Eq. ([TTp. 

In view of the previously discussed relation between the experimental values of r„„ and 

Ta, this equation can be rewritten purely in terms of the latter quantity. In particular, for 

the larger clusters for which Eq. ([13|) holds, we obtain 

1 3r A^^/^ 

When these expressions are used together with d from the previous section in Eq. (^ [note 
that for the smallest clusters it is important to use the full image-charge formula rather than 
the expansion (^], they succeed in describing the behavior of a large fraction of elemental 
metallic clusters. This is shown in Fig. |^ which compares the calculated electron removal 
energies with experimental results on clusters for which data over a sufficient size range are 
available. It is not surprising that the cluster fits are not of uniform qualityill. Indeed, 
whereas our aim has been to reduce as much as possible the number of input parameters, 
the definition of cluster size is inherently not so "universal." The geometric structures of 
clusters containing up to several hundred atoms are neither constant nor necessarily the 



same as the bulk limit (see, e.g., the reviews in Refs. |6^j63|) , and the bond lengths may 
generally vary both with cluster size and from the inner to the outer layers. As a result, the 
formulae given in this section should be viewed as approximations. Their utility is in the 
fact that they provide a sensible interpolation between the atomic and bulk limits discussed 
above, and it is satisfying that a one-parameter expression in many cases provides not only 
qualitative but also quantitative guidance. 

VI. SUMMARY 

The main results of this work can be summarized as follows. 

(1) Starting with an expression for the electron removal energy in terms of the image 
charge potential for an isolated spherical particle and a surface cutoff parameter [Eq. ^, we 
explored how this model may be consistently applied to metallic systems ranging from bulk 
surfaces to finite clusters and down to individual atoms. 

(2) We showed that this approach provides a transparent physical explanation for the 
empirical fact that the atomic ionization potentials and polycrystalline work functions of 
the metallic elements exhibit an almost constant ratio of ~ 2 over the periodic table [Eqs. 

(i 0)]- 

(3) We found that for most elements there is a remarkably close numerical correlation 
between the values of the nearest-neighbor distance in a crystal, the cube root of the atomic 
polarizability, and the image force cutoff parameter. This correlation may be rationalized 
qualitatively, but appears worthy of further study. 

(4) Taking advantage of this correlation, we formulated simple expressions for the cut- 
off distance [Eq. (|^)] and the atomic and cluster radii [Eqs. (|Tl||T^|l^)]. They yield good 



estimates for the work function, the ionization potential, and the cluster electron removal en- 
ergies by using only a single input parameter, the atomic polarizability. No extra adjustable 
parameters are required. 

(5) Generalizing the scaling argument to the case of binary alloys, we found that it can 
simulate the shape of the concentration dependence of work functions of continuous bulk 
solutions [Eq. (P)]. 

(6) We also provided an analytical derivation of the connection between the geometric 
structure of a cluster of close-packed spheres and its outer radius. It is described by an 
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equation of the type ([T3D, but the precise numerical coefficient is shown to depend on the 
packing structure. 
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APPENDIX: HARD SPHERE PACKING 



There exist well-defined geometries for finite systems which provide densely packed struc- 
tures corresponding to specific bulk lattices in the large-size limit. As reviewed in detail, 
e.g., in Refs. p^ , |65| , |66| , these structures are best visualized as arising from the packing of 
spheres. Some of the commonly encountered ones are illustrated in Fig. |^. Here we derive 
the outer radii of these structures as a function of the number of hard spheres, A^. The radius 
is defined as corresponding to the smallest sphere which completely encloses (circumscribes) 
the cluster made up of close-packed hard spheres, each with a radius tq 

1. Cubeoctahedron 

The cubeoctahedron is a shape with small surface area which can be cut out of an fee 
crystal. The number of hard spheres as a function of the number of shells k is given byll 
N = (10/3)A;3 - 5P + (11/3)A; - 1. The outer radius R is given by R = {2k - l)ro. For 
large k this becomes R ~ 2/cro ~ 2(3/10)^/'^roA^^/'^ ~ 1.339roA^^''^. The same expression can 
be derived by evaluating the volume of the circumscribing sphere relative to the sum of the 
volumes of the small hard sphere^!. 

2. Truncated hexagonal bipyramid 

The truncated bipyramid arises from an hep lattice. The number of hard spheres in a 
truncated bi-pyramid is A^ = -(3/4) + (7/2) A; - (21/4)^2 ^ (7/2) A;3 for odd A; and A^ = 
-1 + (7/2) A; - (21/4) + {7/2)k^ for even k. The outer radius is R = {2k - l)ro. For large 
k, R^ 2(2/7)i/3roA^V3 ~ i.3i7roA^i/3. 

3. Rhombic dodecahedron 

The rhombic dodecahedron derives from bcc. The number of hard spheres in a rhombic 
dodecahedron is A^ = Ak^ — 6k'^ -|- 4fc — 1, and the outer radius is i? = (4/ ^/3{k — 1) -|- l)ro 
For large k, R^ 4/y3(l/4)i/3roA^i/3 ^ 1.455roA^^/^ 
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4. Icosahedron 



The icosahedron has the highest symmetry of all discrete point groups. Although due to 
it's five-fold symmetry, the icosahedron does not form bulk crystals, it can be considered 
as a slightly distorted fee crystal. The icosahedron structure has been observed for small 
clusters of inert-gas clusters, Ca, and Mg clusters (see references in ^5|). The relationship 
between the number of hard spheres and the number of shells is the same as for the cube- 
octahedronii: = {10/3)k^ — 5k'^ + (ll/3)/c — 1. The radius is also given by the same 
expression as for the cube-octahedron R = {2k — l)rQ. The resulting expression for the 
cluster radius is therefore the same as for the cube-octahedron. For large k the radius is: 
R ^ 2kro ^ 2{3/10y/^roN^/^ ^ 1.339roA^i/^ 
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Figures 



FIG. 1: A plot of the scaling function ry(d/R) governing the variation of the electron removal energy 
with size, see Eq. 

FIG. 2: Diamonds: the experimental work function^HZl plotted in terms of d = e^/{AW). 
(a)Triangles: r^, the cube root of the atomic polarizability; circles: rnn/'^-, one-half of the crystalline 
nearest neighbor distance, (b) Squares: Eq. ^ for the parameter d. (c) Circles: rQ/2.54. 

FIG. 3: Stars: a plot of the ratio rnn/i'^fa) for the metallic elements (the average is ~ 0.65). 
Diamonds: the ratio IP/WF (the average is ~ 1.8). 

FIG. 4: A plot of metal work functions estimated by Eq. (^) vs. the experimental values. The 
straight lines mark the region of it 10% deviation. 

FIG. 5: A plot of atomic ionization potentials estimated by Eqs. (|8,10) vs. the experimental 
values0. The straight lines mark the region of it 10% deviation. 

FIG. 6: Concentration dependence of the work function of Ag-Au binary alloys. Solid line: estimate 
according to Eq. (ff2|), dotted line: experimental behavior^, dashed line: linear interpolation. 



FIG. 7: A comparison between the experimental electron removal energies of small metal clusters 
(dots) and the model described in the text (solid curves). The first (circled) dot in each plot is the 
atomic ionization potential. The data were adapted from Refs. Q (Li), @|5|,||| (Na), ||4^Ji|] 
(K), II (Y), II (Nb), P (Mn), |2| (Fe), ^] (Co), [|| (Ni), Q (Cu), H (Ag), |6| (Cd), 
(Al), g (In), [|9| (Sn), and [|| (Ce,Pr). 



FIG. 8: Hard-ball and polyhedral depictions of close-packed structures: (a) Cubeoctahedron, (b) 
Rhombic dodecahedron, (c) Truncated hexagonal bipyramid, (d) Icosahedron. 
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